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Ktt S wave scattering is studied using a chiral unitary approach (ChUT) taking 
into account coupled channels. With the amplitudes derived from the lowest order 
chiral Lagrangian as the kernel of a set of coupled channel Bathe-Salpeter equations, 
the I = 1/2 S wave Ktt scattering phase shifts below 1.2 GeV can be fitted by 
one parameter, a subtraction constant, and a scalar resonance corresponding to the 
controversial k (i£g(800)) can be generated dynamically. A good description of the 
I = 3/2 S wave Ktt scattering phase shifts below 1.2 GeV can also be obtained. An 
artificial singularity in the conventional cut-off method of the 2-meson loop integral 
of the ChUT is found. The formalism is applied to deal with the S wave Ktt final 
state interaction (FSI) in the decay J/ip — > K*°(892)K + tt~ , and a qualitatively good 
fit to the data is achieved. The role of k in the decay is discussed. 

PACS numbers: 12.39.Fe, 13.75.Lb, 14.40.Ev, 13.25.Gv 

Keywords: Chiral Perturbation Theory, Bethe-Salpeter equation, Ktt scattering, final state 
interaction 



< Corresponding address. 
*Electronic address: guofk@mail . ihcp . ac .cn| 



2 



I. INTRODUCTION 

The scattering between the lowest pseudoscalar mesons is significant in the under- 
standing of low energy strong interaction and nonperturbative QCD. The most popular 
approach dealing with this problem is Chiral Perturbation Theory (ChPT)Q, 3,0,0,0] 
(for a recent comprehensive review, see The lowest pseudoscalar meson octet 

(tt, K, 77) is identified with the Goldstone bosons associated with the spontaneous 
breakdown of chiral symmetry of QCD. The most general chiral Lagrangian can be 
written in a perturbative manner to an order by expanding the QCD Lagrangian con- 
taining external sources in powers of the external momenta of the pseudoscalar mesons 
and of the light quark masses m u , irid and m s . However, a perturbative expansion to 
any finite order can not describe the appearance of a resonance. For instance, the pure 
ChPT is not suitable for the physics of the scalar isoscalar tttt scattering above 0.6 
GeV where a broad resonance /o(600) (or a) is generally believed to exist. Further- 
more, unitary relation of 5* matrix is only respected by ChPT in a perturbative sense. 
Over the last few years, several nonperturbative methods were proposed to extend 
the chiral expansion to higher energies, such as resummation of loops using a Bethe- 
Salpeter equation (BSE) [2,0|, the inverse amplitudes method (IAM) jil liol llli ll^ 
and dispersion relation methods jLsl. lb"! II^. 

The existence of the scalar particle k (i.e. 7^(800)) and the mass and width of 
this particle, if it exists, has been long controversial since the 1970's HHQ. The 
evidences of the k has been observed in the analysis of the Kit scattering phase shifts 
0, the Dalitz Plot Analysis of the Decay D + -> K~ti + ti + Q, and the BES data of 
the J/i/j decays j^L, 22, 3- 



In order to understand the k, a good understanding of S wave Ktt scattering is 
required. Kit scattering has been investigated using the ChPT by V. Bernard et al. 



at next-to-leading order 24\ and a satisfactory description of the available data up to 
about 900 MeV was found, and by J. Bijnens et al. at next-to-next-to- leading order 25^ . 
The problem has also been investigated using various unitarisation approaches of the 
chiral Lagangian and by other methods such as the interfering Breit- 

Wigner amp litude method [23|, the unitarized meson model and the i^-matrix 
method |29j. A scalar meson corresponding to k was found in all of the mentioned 
unitary chiral approaches. While Long Li et al. found the Kir data could be fitted 
without a pole around 900 MeV in their .fT-matrix approach and only one s-channel 
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resonance between the Kit threshold and 1.6 GeV, i.e., if * (1430) with a mass around 
1438-1486 MeV and a width of about 346 MeV was found 29]. Recently, an important 
work appeared. Some authors have performed a detailed analysis using a set of Roy- 
Steiner (RS) equations j^. The input was high-statistical data at y/s > 1 GeV 
for both 7rK — > 7r.fr and tttt — > KK amplitudes, then the amplitudes below 1 GeV 
were determined. The authors computed S wave scattering lengths and the coupling 
constants Li, L 2 , L 3 and L 4 appearing in the 0(p 4 ) order chiral Lagrangian. The 
precise S wave and P wave phase shifts below about 1 GeV were also obtained by 
solving the RS equations. 

In this paper, we use the on-shell coupled-channel BSE approach (it will be referred 
as ChUT below) proposed by J. A. Oiler and E. Oset :7] to study S wave scattering. 
In Ref. pi, the lowest order ChPT amplitudes are employed as the kernels of a set of 
coupled channel on-shell BSEs. The BSEs are solved to resum the contributions from 
the s-channel loops of the re-scattering between pseudoscalar mesons. The loop integral 
was calculated using a three-momentum cut-off, and this is the only free parameter in 
the approach. By properly choosing the cut-off, the ChUT can well-describe the data 
of the S wave meson-meson interaction up to y/s ~ 1.2 GeV in the isoscalar and 
isovector channels. And more interestingly, the a, /o(980) and ao(980) scalar mesons 
can be generated dynamically as the poles in the second Riemann sheet of the T matrix. 
It was mentioned in Ref. 3J| that satisfactory fits to the Kir scattering phase shifts 



could be produced using a cut-off of 850 MeV and using / = 100 MeV (an average 
between /„. and fx)- However, we find there is a severe problem in the cut-off method 
in calculating the loop integral, which makes it difficult to describe Kn data above 1 
GeV, and we will discuss this in detail in the corresponding section. 

On the other hand, J / ip radiative and hadronic decays provide a rather rich source 
of the information of light hadrons including a and k, and hence a good understanding 
of the interaction between the final state hadrons is required. The ChUT has been used 
to dealing with the meson-meson S wave FSI in 77 — > meson-meson reactions J 3 2j j . in 



some heavy quarkonium decays |33l . |3J, |35j, |36( , and in heavy meson decays 



The paper is organized as follows: The basic definitions and the lowest order chiral 
amplitudes to be used are given in Sec. EJ In Sec. we give a simple review of the 
ChUT first, and then in Sec. IHI Al we discuss an artificial singularity appeared in the 
conventional cut-off method used in the ChUT, and a very simple proof of the unitarity 
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is given in Sec. IIIIBI The calculated S wave Ktt phase shifts are shown in Sec. IIV Al 
and IIV Bl and a good description of the experimental data below 1.2 GeV is achieved. 
The pole position associated to k is given in Sec. IIVD1 As an example of applying the 
ChUT formalism to dealing with the S wave Ktt FSI, the decay J/ip -> K*°(892)K + n- 
is discussed in Sec. El A brief summary is given in Sec. IVII 



II. THE LOWEST ORDER CHPT AMPLITUDES 



A. Definitions 



The lowest order, i.e. 0(p 2 ), Lagrangian for SU(3) ChPT reads 



cm 



(d lt Ud ft rf + M(U + U r )), 



(1) 



where / is the pion decay constant, the physical value of which is 92.4 MeV, and () 
stands for the trace of matrices. The 3x3 special unitary matrix U is defined as 
U = exp and $ is made of the pseudo Goldstone boson fields 



V 



71 

K~ 



71 



K° 



K + \ 
K° 



(2) 



The pseudoscalar meson mass matrix M is given by 

f ml \ 



M 



(3) 



m\ 
v 2m 2 K -m£ J 

where the isospin limit m u = m^ is assumed. The relevant part concerns four mesons 
and can be derived from Eq. (^J) as 

1 



£ (2) 



(4) 



12/ 2 " M 

Let T 7 (s, t, u) denotes the amplitude with total isospin I and 2] (s) its partial wave 
projection with angular momentum I . We take the normalization for partial wave 
amplitudes in such a way that 

1 A 



T!(s) 



dcos9Pi(cos 9)T I {s, t(s, cos 9), u(s, cos 9)), 



(5) 
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where Pi (cos 9) is the Z-th order Legendre polynomial, 9 is the scattering angle in the 
center of mass frame. The Mandelstam variables s — (pi + p 2 ) 2 , t — (p± — p^) 2 and 
u = (pi — Pi) 2 are bound to the identity s + t + u = Ylt=i m \ when the particles are 
on-shell, and 

t(s, cos 9) = ml + ml - 2^[m\ + X ^ m "\ ][mj ~ ^g^g ] 

+ —\j X(s, ml, m|) X(s, m|, m|) cos 9 (6) 

with A(s, m 2 , m 2 ) — [s — (mi + mj) 2 ][s — (mj — mj) 2 ]. 

We shall use the phase convention \ir + ) = — 11,1) and |^°) = — 1 1/2, 1/2). The 
calculations of / = 1/2 scattering will be performed with single channel (Kir), two 
channels (Kit and Krj) and three channels (Kir, Kr\ and Kr]'), respectively, and the 
single channel 1 = 3/2 Kir scattering will also be investigated. The lowest order chiral 
amplitudes to be used in our coupled channel calculations of the S wave Kn scattering 
will be given in the following subsections. We will label them with V and label the full 
amplitudes in the BSEs with T to avoid confusion. Since we will restrict ourselves in 
the S wave, the subscript I will be neglected in the following and I = will be implied 
everywhere. 

B. Two channels case 

In this case, we take rj = r]g and the mixing between r\ and rj will not be considered. 
Let us label the channels and Kt] with 1 and 2 respectively, such that t,u), 
Vi 2 (s, t,u) and V^s^t^u) represent the amplitudes with total isospin / for Krj — > Km 
Krj — > Kn and Krj — > Km By time reversal invariance, one has V/j(s, t, u) = V^(s, t, u). 

ir + K + — > tt + K + is a pure 1 = 3/2 process, 

V^ 2 (s,t,u) = V w+K+ (s,t,u). (7) 

1 /2 

According to the crossing symmetry, the isospin relation for V x { (s, t, u) and 
Vi2 2 (s, t, u) can be interpreted as 

V^(s,t,u) = \v^ 2 (u^s)- l -V'/ 2 {s^u\ (8) 
Vl 2 /2 (s,t,u) = V3V K o ri ^ K o 7r o(s,t,u). (9) 
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Now the 0(p 2 ) order ChPT amplitudes can be derived from Eq. (@J). One has a pure 
S wave amplitude for 7 = 3/2 process 

1 , 



V 



2K s-ml-m 2 K ), 



and for I = 1/2 the amplitudes contain both S wave and P wave 



(10) 



Vl/ 2 (s,t,u) 



4/ 2 

1 



(4s + 3t - 4m 
1 



4m 



K)i 



3t+ -ml + -m 2 K + m 2 v ), 



4/ 2 

1 / 2 2 
r(3*-;X 



2mf 



(11) 
(12) 

(13) 



4/2-" 3 

Note that the calculations will be performed using on-shell BSE, so that the on-shell 
condition p 2 = m 2 has been used in writing the above amplitudes. The isospin phase 
convention of K is different from the one in Ref. Q, such that there is an extra 
minus before Eq. (I12H compared with T2(s,t,u) in Eq. (B14) therein, but the physical 
observables are not influenced by different phase conventions as they should be. 



C. Three channels case 

If one want to see the possible impacts of the Krj' channel, r/ should be included in 
the chiral Lagrangian. One way is extending the SU(3) matrix to a U (3) one 



nonet 



* + 7=3 V0 



V2' 



\ 



71 

K~ 



where the standard rj — rf mixing 



TV 



+ 



V3 

K° 



\v) = 3^0) + 



2V2 



W) 



2V2 



\Vo) ~ 77 \Vs) 



k+ ^ 

K° 



(14) 



(15) 



is assumed. We label the channel Krj' with 3, such that V^ 2 (s,t,u), V^' e (s, t, u) and 

1 /2 

V 3 s (s,t,u) represent the amplitudes for Krf — > Kir, Krj' — > Kr\ and Kr( — > Krj', 



rl/2. 
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respectively. The amplitudes concerning only K and n do not change while the others 
change to 



V&\8,t,u) = -j^(-3t + 2m 2 K + m 2 v ), (16) 
V 2 T(s, t, u) = -gj5 (3t -m 2 K - 2ml). ( 17 ) 



The amplitudes involving rf are 



vU\s, t, u) = ^(-3t + 3ml + 8m 2 K + mj), (18) 



1 

12/2 

V^ 2 (s, t, u) = ^(3t - 3m 2 + 2m 2 , - m 2 - m 2 ,), (19) 
7i /2 («, f, u) = (* - 2m 2 + 12m 2 , + |m»,). (20) 

III. COUPLED-CHANNEL CHIRAL UNITARY APPROACH 

In our normalization, the unitary relation for the partial amplitudes with isospin / 
satisfies 

ImT\s) = -T I \s)p(s)T I (s), (21) 

where T 7 (s) and p(s) are nxn matrices for n channel calculations. The matrix elements 
T*Js) will be given by the on-shell BSEs. p(s) is a diagonal matrix with 

*M - (22) 

where p cmi is the three-momentum of one meson in the center-of-mass (cm.) frame for 
the i-th channel. This momentum is given by 

Pcm = T^yV - (mi + m 2 ) 2 ] [s - (mi - m 2 ) 2 ] (23) 

in a channel with two mesons of masses m x and m 2 respectively. In what follows we 
will omit the label I for simplication. 

Following Ref. jlj], the amplitude of meson-meson scattering can be cast using the 
BSE 

T = (1 - VG)~ X V, (24) 

where V is the lowest order chiral amplitude and G is the loop propagator. For coupled 
channel calculations, all of T, V and G are matrices. One important feature of this 
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approach is that V can be factorized on shell in the BSEs, and so that the integral 
equations become algebraic equations which can be solved simply. This feature is 
justified using the N/D method of dispersion relations and by a comprehensive 
treatment of the BSE 



A. Artificial singularity in the cut-off method 

The loop integral in the z-th channel is 

r d 4 q 1 1 

ll[s) = % J (2tt) 4 g 2 - m\ + is ( Pl + p 2 - qf - m\ + is' (25) 

where p\ and p 2 are the four-momentum of the two initial particles respectively, rri\ 

and m2 are the masses of the two particles appearing in the loop. Sometimes in ChUT 

the integral is calculated by using a three-momentum cut-off parameter q max in the 

cm. frame |7[. The analytic expression can be worked out as 

/-f/ % 1 r / S + A S ~ A ^ 

G (s = 7— icriarctan — h arctan — - — 

_ [(s + A) ln(^(l + AO) + (s — A) ln(^(l + A 2 ))]} (26) 

777i 7772 



where a = [— (s — (mi + 777 2 ) 2 )(s — (mi — t^) 2 )] 1 / 2 , A = m 2 — mj and Aj = 
(i = 1, 2), and the channel label ii has been dropped. 

However, there are some problems. If one goes to values of the on-shell momenta 
above the cut-off, one can not use the cut-off formula j^. On the other hand, we find 
that there is a severe problem using this cut-off method to calculate the loop integral. 
From Eq. (|26|). one can find there is a singularity above the threshold which is located 
at 



(Vising = V m l + l^nax + V m 2 + llax' ( 27 ) 

When y/s goes to (\fs) singi the real part of G(s) will goes to infinity, and when yfs 
is greater than (\fs) sing given below in Eq. (|27)l . Eq. (|2^j) gives a vanishing imaginary 
part. However, if momenta are large enough to bring the propagators on shell, the 
imaginary part should be given by — p cm /(87r-y/i) . We can give some remarks on the 
singularity. First, after some manipulations, the real part of the loop integral can be 
given as 

ReG W = >P r- i\ q \ f""+^ (28) 



9 



where P / represent principal integral, and cjj = \Jq 2 + m 2 . One can see that the 
singularity comes from the possibility of the denominator in the integrant being zero. 
And when the cut-off q max runs to infinity, the singularity will disappear. Second, the 
singularity rises from the arctan/(s) function. If we perform a series expansion at 
the singularity point, the first term should be arctanz, or equivalently ^ InO using the 
identity arctan z = ^ In ^=4. It should be a logarithmic singularity. Third, when y/s = 
(y/s)sing, the three-momenta of the propagators equal to q ma x- Then no momentum 
is available to bring y/s larger than (y/s) sin g, and the imaginary part of G(s) should 
vanish as it appears in Eq. (|2fijl . In other words, (y/s) S i ng should be available limit 
of this method. On the other hand, there is still a region below (y/s) S i ng where the 
impact of the singularity is large, and one should be care of that. For the S wave Kit 
scattering, a typical value of the q m ax is about 0.6 GeV from fitting to the data when 
the kernel of the BSE is of the 0(p 2 ) order. The corresponding position of the artificial 
pole will be about 1.4 GeV, and its impact will show up above 1 GeV apparently. Then 
it is difficult to describe Kit data above about 1 GeV in this formalism. Note that this 
cut-off method is different from the one used in Pauli-Villars regularization (see e.g. 
Ref. 4oT|). In this method, the cut-off is just a three-momentum cut-off. 

The way to solve this problem is using a dimensional regularization method with 
a dispersion relation to deal with the loop integral. The analytic expression has been 



given as 



41] 



2 



„, N 1 , , . , mi A — s . mi 

GW =ifa5 {a( " )+lof v + — log 4 

+ |^[log (s - A + (x) + log (s + A + a) 

-log(-s + A + o-) -log(-s- A + o-)]}, (29) 

where a(fi) is a subtraction constant, fi is the regularization scale. The result is inde- 
pendent of \i because the change causing by a change of \i can be cancelled by a change 
of the subtraction constant a(/x). 

In order to show the impact of the artificial singularity when the cut-off method is 
used, we plot the real part of the loop integral with m\ = m n and m2 = mx in Fig. ^ 
In Fig.HJ the solid curve represents the result from the cut-off method with q max = 0.57 
GeV and the dashed curve represents the result from the dimensional regularization 
method calculated at fi = mx with a(rrik) = —1.41. From Eq. (|27|) . the singularity is 
located at yfs = 1.34 GeV, and the values above 1 GeV are all affected significantly. 
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FIG. 1: The real part of the loop integral calculated using the cut-off method (solid surve) 
and the dimensional regularization method (dashed curve). 



However, the physical results of Ref. |7[ where the cut-off method was used will not 
be influenced, because with the parameter q max = 1.03 GeV therein the singularity 
will be located 2.08 GeV, 2.29 GeV and 2.21 GeV for the vrvr, KK and irr] channels 
respectively, and the effect below 1.2 GeV is very small. 



B. A simple proof of unitarity 

It has been shown that by this coupled-channel BSE approach unitarity is ensured 
jloL lol ]. In fact, it can be realized quickly as follows: 
The unitary relation Eq. ()21j) can be rewritten as 

1 1 



T(s)* T(s) 
The general solution of Eq. (|30|) is of the form 



2ip(s). (30) 



K(s) + ip(s) 

where K(s) is a real function of s. While by the Cutcosky rule j^, the imaginary part 
of the loop integral can be related to the phase space factor as 

ImG(s) = -p(s). (32) 

So Eq. (|24j) can be rewritten as 
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1 

~ V~ 1 (s) - ReG(s) - ilmG(s) 

= K(s) + i P ( s y (33) 

where K(s) = V _1 (s) — ReG(s). In this approach, V(s) is the lowest order ChPT 
amplitude and so that is real. From the above equation one can see in this approach 
Eq. ({23| ensures unitarity actually. Furthermore, if one use such a kernel that the BSE 
can be dealt with in an on-shell way, the unitary relation can be ensured if and only if 
the kernel is real. 



IV. RESULTS AND DISCUSSION 

A. S wave 1 = 1/2 Ktt phase shifts 

Given the above knowledge, one can calculate the S wave Ktt scattering phase shifts 
numerically. The only parameter in our calculations is the subtraction constant a([x). 
In principle, it should be different for different channels and be same only under exact 
SU(3) symmetry. In practice, in order to reduce the number of free parameters, we 
take a(/i) to be same for different channels. All of the results presented in this paper 
are calculated at [A = rriK- In order to get the most accurate description, we use the 
MINUIT function minimization and error analysis package from the CERN Program 
Library j^J to find the most appropriate value of the free parameter ajrriK 



The S wave 7 = 1/2 Kit phase shift data are taken from Refs. 



4J,|45|, |46j. The 



parameters can be determined from fitting to the data, and the values are 

for one channel : a(m^) = —1.412 ± 0.017, 
for two channels : a(mj:) = —1.278 ± 0.014, 
for three channels : ct(mx) = —1.383 ± 0.006. (34) 

The calculated phase shifts are plotted as solid curves in Fig. El where the dashed line 
depicts one-channel case, the dotted line depicts two-channel case, and the solid line 
depicts three-channel case, respectively. One can see the data below 1.2 GeV can be 
well described for one channel and two channels, and the data below 1.2 GeV can be 
well described for three channels. The fit from three-channel calculations is slightly 
better. 
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FIG. 2: S wave 1=1/2 Ktt phase shifts with the dashed line depicting one-channel case, 
the dotted line depicting two-channel case, and the solid line depicting three-channel case, 
respectively. 



As mentioned in the previous section, if one use the three-momentum cut-off method 
to deal with the loop integral, there will be an artificial singularity in the loop function. 
The typical value of q m ax is 0.4-0.6 GeV from fitting the phase shift data below 1 GeV 
or 1.2 GeV, and the singularity is located at about 1.0-1.4 GeV. Then the effect of the 
singularity above 1 GeV is large, see Fig.HJ That is to say, one cannot get a satisfactory 
description of Kir S wave scattering above 1 GeV using the three-momentum cut-off 
method. One way out is using the dimensional regularization to deal with the loop 
integral as illustrated above. 



B. S wave 1 = 3/2 Ktt phase shifts 



The experimental data of S wave / = 3/2 Ktt phase shifts are taken from 



Refs. 



1«. 



46]. The calculated results are plotted in Fig. El From fitting to the data, 
the subtraction constant is determined to be almx) = —4.643 ± 0.083. 



C. Comparison with the results from Roy-Steiner equation analysis 



Very recently, some authors have performed a complicated analysis of Ktt scattering 
using RS equations j^. The experimental data above s = 0.935 GeV 2 were used as 
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FIG. 3: S wave / = 3/2 Kir phase shifts. 



input, and the solutions were obtained in the range s < 0.935 GeV 2 . The results are 
very precise, and it is deserved for our results to be compared with the RS equation 
results. For illustration, we plot our results in the range s < 1 GeV 2 for 1=1/2 with 
two channels and for / = 3/2 as well as the RS equation results [3^ in Fig. 0] The 
upper half and the lower half of Fig. H]are for 1=1/2 and / = 3/2, respectively. In 
this figure, the solid lines represent the ChUA results, and the dashed lines represent 
the results from the detailed analysis using the RS equations j^jj]. From comparison, 
we find that our results are consistent with the precise results in Ref. [3^] below 1 GeV 
for both of I = 1/2 and I = 3/2. 



D. Poles in the amplitudes 

The physical resonances can be associated to the poles in unphysical Riemann sheets 
of the scattering amplitudes. After the same analytic continuation of the loop function 
G(s) as in Refs. 0, 0|> one can find a pole corresponding to the long controversial 
scalar particle k which is called 7^(800) in the PDG [48]. The pole position for each 
calculation is given in Table |1J 

In order to see the stability of the pole position against different regularization 
schemes, the results from the cut-off method are also listed. For each calculation, 
the parameter q max is determined from fitting to the phase shift data below 1 GeV 
where the influence of the artificial singularity is small as pointed out in the previous 
section. In the first column of Table HI DRM represents the dimensional regularization 
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FIG. 4: Comparison of S wave Kir phase shifts with the results from Roy-Steiner equation 
analysis. The upper half and the lower half are for / = 1/2 and / = 3/2, respectively. The 
solid and dashed lines depict our results and the RS equation results, respectively. 



method and COM represents the cut-off method. In the second column of Table U the 
parameter is a(mx) for DRM and is q m ax for COM. The pole positions are calculated 
using the central value of each parameter. One can see that the difference of pole 
position, especially the real part, is small using different regularization schemes. 

TABLE I: Pole positions found in the second Riemann sheet of the T matrices. 



Method Number 


of channels Parameter 


Pole position (GeV) 




1 


-1.412 ±0.017 


0.736 - 


i0.273 


DRM 


2 


-1.278 ±0.014 


0.725 - 


i0.297 




3 


-1.383 ±0.006 


0.742 - 


i0.273 




1 


(0.512 ± 0.005) GeV 


0.729 - 


i0.244 


COM 


2 


(0.477 ± 0.003) GeV 


0.730 - 


i0.244 




3 


0.486 GeV 


0.728 - 


i0.247 



Because in the lowest order, the decay constants are same for tt, K and r/, one 
can investigate some effects of higher order Lagrangian on the pole position by taking 
the physical values of the decay constants of K and i], that is, fx = 113.0 MeV and 
f v = 110.9 MeV (4^. Let us study the case of two channels, and the results of the other 
cases are similar. The replacement can be done following the Appendix B of Ref. Tp| . 
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For example, 1//^ in Eq. (fTTJ) will be replaced with The pole positions are 

listed in Table |TI] The real and imaginary parts of the pole position are both slightly 
larger than the corresponding ones listed in Table HI 

TABLE II: Pole positions found in the second Riemann sheet of the two channel T matrices 
taking different decay constants for it, K and rj. 

Method Parameter Pole position (GeV) 
DRM -1.798 ±0.015 0.790 - i0.320 
COM 0.619 ± 0.006 GeV 0.773 - i0.276 



These results can be compared with previous results in other chiral unitary ap- 
proaches. The pole position appeared in the second Riemann sheet associated to k 
were found to be around 770 — z250 MeV in Ref. 0] where COM was used and 
779 + i330 MeV in Ref. where DRM was used. These results are obtained using 
the 0(p A ) chiral amplitudes as the four- meson contact terms. The position of the pole 
we found is similar to the one found in other approaches. Furthermore, from com- 
parison, it seems that the higher order corrections from the decay constants give an 
important contribution to the 0(p 4 ) order corrections to the pole position. 



V. APPLICATION TO THE DECAY J/ip -» i^* (892)K + vr 

In the preceding sections, we have constructed a formalism of the S wave Kit scat- 
tering below 1.2 GeV. In this section, as an example of its applications, it will be 
applied to describe the S wave Kn FSI in the decay J/ip — > K* (892)K + n~ . We will 
calculate the invariant mass of K + n~ below 1.2 GeV and fit to the experimental data 
The relevant decay mechanisms are plotted in Fig. 

In Fig. (a) represents the direct decay term, and (b) takes into account the S 
wave K7r FSI with coupled channels. We are interested in the region below 1.2 GeV. 
In this region, only one resonance, K*°(892), can decay into K + n~ . However, the 
branching ratio of J/ip — > i^*°(892)^*°(892) is very small. As a simple consideration, 
no intermediate resonances will be considered. 
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K'/ J' 



J/W K*° J/W K*° 

(a) (b) 



FIG. 5: Diagrams for the decay J/ip A*°(892)A+7r 



Similar to Refs. the contact terms of the vertex J / ipVPP can be extracted 

from a Lagrangian of the form 

£ c = g^(V^$), (35) 

where ip^, and $ are the fields of J/ip, vector and pseudoscalar octets respectively. 
V and $ are both 3x3 matrices with 



V 



\ 



K* 



AT* / 



(36) 



where ideal mixing between and u is assumed, and $ has been given in Eq. (j2J. The 
relevant terms in Eq. (J35j) can be written as 

(37) 

Let us denote the momenta of J /if), K + , n~ and A'* (892) with p, p 1 , p 2 and p 3 , 
respectively. From Eq. (|3*7j) . the amplitude for Fig. Efa) is of S wave, and it is 

t = -ge£\p)e^'\p 3 ), (38) 

where el^(p) and e^ x '\p 3 ) are the polarization vectors of J /if) and A'* . 

From the preceding sections, the S wave 7 = 1/2 AV phase shift data below 1.2 GeV 
can be fitted well using two coupled channels. The coupled channel S wave Kir FSI can 
be taken into account using the formalism. After that, the amplitude corresponding 
to Fig. [31(a) and (b) is 

t c = -ge^{p)e^{p 3 )F{s 12 ) 



ge M(p) E vWfa)(l + Gll (K + Tr-\t\K + Tr-) 



V2 



Gn(K + 7r-\t\K°ir ) 



V6 



G 22 (K + n-\t\K°r ] )), 



(39) 
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where (K + K~\t\P\P2) (-P1-P2 = K + tt~ , etc.) represents the full amplitude of the process 
P\Pi — > if + 7r~ calculated in the ChUT, and S12 = (pi + P2) 2 = tu^. Because the 
isospin of J /if) is and of K*° is 1/2, the FSI amplitudes must be of isospin 1/2. These 
amplitudes can be related to the amplitudes in the isospin basis as 

(K + *-\t\K+*-) I=1 ft = ? i T}p(8), (40) 
(K + n~\t\K°n°) I=1 / 2 = -^^(a), (41) 



\s), (42) 



where T^ 2 \s) is the full amplitude of the channel j — *■ i (i — 1 : if-7r, z = 2 : if 77). 
The invariant mass spectrum of K + n~ can be calculated using the formula [48] 

dT 1 1 F ££/l*clVillPs|<ffW, (43) 



dmi 2 (2tt) 5 16M 2 



A A' 



where X)a Sa' describes the average over initial states and the sum over final states, 
(|Pil,^i) is the momentum of K + in the cm. frame of K + tt~ , and (Ipal,^) is the 
momentum of _K"*°(892) in the rest frame of J/if>. 

On the other hand, because the experimental data are without sideband subtraction 
we will consider the sequential process J /if} — * K*°(892)K~tt + — > K + n~ K~n + as 
a sideband background. 

For the process J /if} — > K*°K~ir + , one can write an effective amplitude like the one 
in Eq. Ij59}l . 



geW(p)e^(p*°)F(s u ), (44) 



where s 34 = (p K - +p 7r +) 2 = m\±. 

The amplitude of the P wave decay if* (892) — ► K + n~ can be written as 

t K *o K + w - = g K *K^ x "\p A ) ■ (p 2 - pi). (45) 

From the width of if* (892), the dimensionless coupling constant can be determined 
to be gK*KiT — 4.61. Now one can write the amplitude for the sideband background 
process as 

rp, \ (AW ^ U + {P1+P2Y{PI+P2) U /S 12 ( 

t s = -ggK*K^F{s^)e^\p) — , — — - — — (p 2 - pi)„, (46) 
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where Tk*{s\2) is the momentum-dependent width of K*°(892). The relation of 
^k*(su) and the on-shell width Tk*(itVk*) = 50.7 MeV is 

T K ,(s 12 )=T K 4m 2 K ,)^^, (47) 

S12 go 

where g and go is momentum of one meson in the cm. frame 

q = 7 — v/ (si2 - K + ^) 2 )(gi2 - (m K - m n ) 2 ), (48) 



g = 77^ — J {m 2 K , - (m K + m n ) 2 )(m 2 K , - (m K - m w ) 2 ). (49) 

Since K + and ir~ in their cm. frame are in S wave in the direct decay process and 
in P wave in the sideband background process, there is no interference between the 
amplitudes t c and tb- Hence the invariant mass of K + ir~ can be calculated using the 
formula 3] 

dT 1 1 



dm K * (2tt) 5 16M 2 



EE / |t c | 2 |p!i|p 3 |^^3 



A A 



+ j^ys^p j J a dm M ^\tf\p{\\p';\dnidn' 3 *dn 12 (so) 

where J2x Sa' describes the average over initial states and the sum over final states, 
(|Pi|,^*) is the momentum of K + in the cm. frame of K + n~ , (|p 3 |,f2 3 ) is the mo- 
mentum of ^*°(892) in the rest frame of J/if), and (|p' 3 *|, ^3*) is the momentum of 
K~ in the cm. frame of K~tt + . In the integral J a dmsA in Eq. fl5()(l . a = 0.812 GeV 
and b = 0.972 GeV according to the experimental events selection criterion 0.812 GeV 
< m K - w+ < 0.972 GeV U. 

n 

There is only one parameter g. The program package MINUIT |42J is used to make 
a fit and find the most appropriate value of the parameter. The BES data are taken 
from Ref. j^j. The results are plotted in Fig. El and the fit is good qualitatively. 
The solid line represents the fit result, the dash-dotted and dashed lines represent the 
contributions of the direct term and the direct term + the S wave Kn FSI, respec- 
tively. The dotted line represents the background or sideband contribution coming 
from J/ij -> K*°K-n + -> K + it-K-it + in the region 0.812 GeV < m K - n+ < 0.972 
GeV. The dash-dot-dotted line represents the contributions concerning k or of the FSI 
and the interference of the FSI and the direct term. One can see the wide bump below 
1.1 GeV besides the peak of K*°(892) is due to k. 
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-200: 



0.6 0.7 0.8 0.9 1.0 1.1 1.2 
M KV (GeV) 

FIG. 6: The invariant mass spectrum of K + tt~ below 1.2 GeV of the decay J/tp — > 
K*°(892)K + ir~ . The black dots with error bars are the experimental data. The solid line 
represents the fit result, the dash-dotted, dashed and dotted lines represent the contributions 
of the direct term, the direct term + the S wave Ktt FSI and the sideband background process 
J/ip — > K*°(892)K~tt + — > K + tt~ K~tt + , respectively. The dash-dot-dotted line represents 
the total contributions of the FSI and the interference of the FSI and the direct term. 

A decay amplitude can be expressed as a sum of a background and some Breit- 
Wigner amplitudes associated to the relevant resonances. From the experimental data 
plotted in Fig. El one can see that there is a wide bump below 1.1 GeV besides the peak 
of if* (892). Were there no k and were the background of no structure and just like 
a form of Eq. ()38)1 . this pump would be hardly understood. From our results plotted 
in Fig. IHl the wide bump is generated by the S wave Ktt FSI. Note that in our ChUT 
approach, the scalar particle k is generated dynamically by the Ktt FSI. That is to see, 
the bump is generated by k. 

The value of the parameter determined from fitting to the spectrum of the number 
of events is g = 1770.3. However, it is not physical, and so far we cannot determine 
its physical value because the branching ratio B(J/ip — > K*°(892)K + tt~~) has not been 
known yet. 

VI. CONCLUSIONS 

We have used a chiral unitary approach to study S wave Ktt scattering and gotten 
a satisfactory description of phase shift data below 1.2 GeV. On-shell BSE is used, and 
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the lowest order ChPT amplitudes are used as its kernel. One-channel, two-channel 
and three-channel calculations are performed for the 1 = 1/2 Ktt scattering, and one- 
channel for / = 3/2. There is only one parameter a(mjf) for each isospin in the model, 
and its value is determined from fitting to the data. A pole of the I = 1/2 T matrix 
is found in the second Riemann sheet. We investigate the pole position with different 
number of channels with both cut-off and dimensional regularization schemes, and find 
that the pole position is stable against different regularization schemes. If we take into 
account some higher order corrections, the real and imaginary parts of the pole are both 
slightly enhanced. This pole corresponds to the long controversial broad resonance k. 
In this way, the mass and width of n are estimated to be about 757 ± 33 MeV and 
558 ± 82 MeV. It is pointed out that there would be a singularity in the loop function 
if the conventional cut-off method were used. The singularity comes from that the 
three-momentum cut-off set a limit for the applicable region. 

One advantage of this approach is that there is only one parameter. Once this 
parameter is determined from the phase shift data, it can be used to other processes 
such as heavy particle decays involving S wave Ktt or Kt] FSI. In recent years, there 
are other excellent works on Ktt scattering, such as In Ref. [2q , the resonance 

chiral Lagrangian together with a unitarization was used, and the data of S wave 
Ktt scattering below 2 GeV can be fitted well. However, the number of free parameters 
is large. The pole position corresponding to the k is about 700 — «300 MeV. It is similar 
to the one found in our work. The difference of the mass is large slightly. However, 
the width of the resonance is about 600 MeV, and this makes a small difference of the 
m ass and width is of little hnportanee. In Ref. 0, the anthers solved a set of six 
RS equations for the S and P waves of the Ktt scattering amplitudes. The input was 
high-statistical data at ^/s > 1 GeV for both ttK — > ttK and tttt — > KK amplitudes, 
then the amplitudes near threshold were determined. Certainly, this is an important 
work, and we compared our results for S wave Ktt scattering phase shifts below 1 GeV 
with the results in this work. From comparison, we find that our results are consistent 
with the precise results in Ref. jQ. Furthermore, the one-parameter model can be 
used easily to deal with S wave Ktt FSI appeared in other processes below iurk — 1-2 
GeV. As an example, we calculated the invariant mass spectrum of K + tt~ of the decay 
J ftp — > K*° '(892) K + tt~ , and the fit to the experimental data is qualitatively good 
below 1.1 GeV. From the analysis of different contributions from different processes, 
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we found that k is important to produce the wide bump below 1.1 GeV. 
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